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ABSTRACT 

Equations of fully general relativistic radiation hydrodynamics in Kerr space-time 
are derived. While the interactions between matter and radiation are introduced in 
the comoving frame, the derivatives used when describing the global evolutions of 
both the matter and the radiation are given in the Boyer-Lindquist frame (BLF) 
which is a frame fixed to the coordinate describing the central black hole. Around a 
rotating black hole, both the matter and the radiation are influenced by the frame- 
dragging effects due to the black hole's rotation. As a fixed frame, we use the locally 
non-rotating reference frame (LNRF) which is one of the orthonormal frame. While 
the special relativistic effects such as beaming effects are introduced by the Lorentz 
transformation connecting the comoving frame and the LNRF, the general relativistic 
effects such as frame-dragging and gravitational redshift are introduced by the tetrads 
connecting the LNRF and the BLF. 

Key words: accretion: accretion disks — black hole physics — hydrodynamics — ra- 
diative transfer — relativity. 



1 INTRODUCTION 

Accretion disks around black holes play one of the essential roles in the active universe. When the mass accretion rate is 
near or over the Eddington mass accretion rate, the radiation in the disk interact with matter. For the supercritical accretion 
flows considered in e.g. Syfert galaxies or the black hole binaries, the photons interact with matter. On the other hand, for 
the hypercritical accretion flows considered as one of the leading candidate for the central engine of gamma-ray bursts, the 
neutrinos interact with matter around the central black hole. In such situations, the dynamics and the energy balance in 
matter and radiation are influenced by each other. The radiative transfer in the curved space is investigated by many authors 
(e.g. Lindquist 1966; Anderson & Spiegel 1972; Schmid-Burgk 1978; Thorne 1981; Schinder, B ludman & Piran 1988; Turolla & 
Nobili 1988; Anile & Romano 1992; Cardall & Mezzacappa 2003; Park 2006b iLindauistJ [| 19661 ) gives a general treatment of the 
radiation transfer equation and the radiation hydrodynamic equations derived from the zeroth and the first moments of the 
radiat ion t ransfer equatio ns by using a comoving Lagrangian fra me of refe r ence. For the higher moments, lAnderson fc Spiegel 
1 1972h and Thorne 1 1981 ) gives a detailed treatments. Especially, Thornel 1 198lh derive general relativistic moment equations 
up to an arbitrary order by introducing projected symmetric trace-free (PSTF) tensors. This formalism has been used in 
spherically symm etric problems ( e.g. Flammang 1982, 1984: Turol la & Nobili 1988). Other forms of the moment equations 
are presented by ISchinder (|l988h and ISchinder fe Bludmanl (|l989h in a static, spherical space-time by using a Lagrangian 
comoving coordinate system. In these formalisms based on the comoving frames, the physical quantities in terms of matter 
and radiation and the di rectional deriva tives are described in the comoving frame. 

On the other hand, Mihalasl 1 1980[ ) introduce the alternative approach by using the radiation hydrodynamic equations 
in the Eulerian framework where the derivatives become much simpler form in this framework. In this approach, the physical 
quantities in terms of matter and radiation and the derivatives were introduced in the frame fixed to the coordinate of the 
central object, e.g. a black hole, while the interactions between matter and radiation were calculated in the comoving frame. 
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That is, the local processes like the interaction between the matter and the radiations are evaluated in the comoving frame, 
while the derivatives which are used for the calculations of the global dynamics of both the matter and the radiation are 
derived in the frame fixed to the coordinate describing the central object. In principle, t his fo r malis m can be extended to 
arbitrary space-time. Actually, by using the energy- momentum tensor i n a covarian t form, Par id 1 19931 ) derived the radiation 
hydrodynamic equations for spherically symmetric systems, and recently, Park ( 20061 ) give the explicit expressions for the basic 
equations of the general relativistic radiation hydrodynamics in Schwarzschild space-time. So far, the formalism in the space- 
time around a rotating black hole which is usually considered for the black holes in the real world have not been presented. 
One of the natural next step is to derive the explicit expressions for the general relativistic radiation hydrodynamics in space- 
time around a rotating black hole in this formalism. In this paper, we assume c = 1 in most equations except in a few cases 
where c is explicitly used for clarity. After giving the general forms of the basic equations for the general relativistic radiation 
hydrodynamics in §2, in §3 by using the orthonormal tetrads fixed to the coordinate (§3.1), the radiation moments (§3.2), 
the radiation four-force (§3.3) and the radiation hydrodynamic equations (§3.4) in Kerr space-time are derived. Concluding 
remarks are given in the last section. 



2 BASIC EQUATIONS FOR GENERAL RELATIVISTIC RADIATION HYDRODYNAMIC IN 
COVARIANT FORM 

Here, we summarize the basic equations for the general relativistic radiation hydrodynamic in covariant form (e.g. Mihalas & 
Miharas 1984). 

The energy-momentum tensor for matter of an ideal gas is described as 

T<* 9 = p o h g u a u + P e g a(3 , (1) 

where u a , po, h s and P g are the four- velocity, the rest-mass density, the relativistic specific enthalpy and the pressure of the 
gas, respectively. The specific enthalpy of the gas, h s , is calculated as h g = (e g + P s )/po where e g is the energy density of 
the gas. Here, the fluid quantities h g , e g , P g and po are all being measured in the comoving frame of the fluid. The radiation 
stress-energy tensor is given as 

^ = / / 7< ' n ' V ^ na ^ idvdQ -> ( 2 ) 

where I(x a ; n, v) is the specific intensity of photons moving in direction n with the frequency v. The photon four- momentum 
p a = (hv)(l,a) and n a = p a /hu. 

The continuity equation, i.e. particle number conservation equation, in the absence of particle creation and annihilation 
is given 

{nu a ). a = 0, (3) 

where n is the particle number density measured in the comoving frame. Here, n is related to po as n — po/m g where m g is the 
mass of the gas particle. On the other hand, the conservation equation for the total energy-momentum of gas plus radiation 
is given as 

(T aS + R a9 ). = 0. (4) 
The radiation four-force density acting on the matter is given as 

G a — i J J [xvl(n, v) — rj u ]n a 'dvdfl, (5) 

where \v an< i 1v are the opacity and the emissivity, respectively. The invariant emissivity and invariant opacity are r/^/u 2 
and vxv , respectively. The dynamical equations for the matter and the radiation field are described as 

T Q/3 ;/3 = G a , (6) 
and 

R af> ifi = _ G c. i (7) 
respectively. 



3 RADIATION HYDRODYNAMICS IN KERR SPACE-TIME 

In order to see the correspondence with the work by Park ( 20061 ) deriving the equations for the radiation hydrodynamics in 
the Schwarzschild space-time, in the present study we use the Boyer-Lindquist coordinate for the description of the rotating 
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black hole's space-time. The background geometry around the rotating black hole written by the Boyer-Lindquist coordinate 
is described as 



ds 2 = g ol /3dx a dx 13 , 



■ -ii j {dx i + i dt){dx i + /3 j dt), 



(8) 



where i, j = r, 9, <fi and the nonzero components of the lapse function a, the shift vector (3 l and the spatial matrix 7^ are 
given in the geometric unit as 



Asm 2 9 



A' 1M 



14>4> 



(9) 



Here, we use the geometric mass m = GM/c 2 



+ a 2 cos 2 9, A = r 2 - 2Mr + a 2 and A = (r 2 + a 2 ) 2 - a 2 A sin 2 9 = 



EA + 2mr(r 2 + a 2 ), where M is the black hole mass, G is the gravitational constant and c is the speed of light. The position 
of the outer and inner horizon, r±, is calculated from A = as r± = m ± (m 2 — a 2 ) 1 ^ 2 . The angular velocity of the frame 
dragging due to the black hole's rotation is calculated asu = — flt^/ ga^ = 2ma r/A. When a/m = 0, we obtain a = F, {3r — 0, 
7 rr = 1/r 2 , ■yee = r 2 and 7^ = r 2 sin 2 9 where fs(l- 2m/r) 1/2 [|Parkll2006l ). 



3.1 Orthonormal Tetrads 



In the present study, we use three frames: (1) the Boyer-Lindquist frame (BLF) which is the frame based on the Boyer- 
Lindquist coordinate describing the metric and have been frequently used for the description of the global dynamics (e.g. 
Gammie & Popham 1998), (2) the locally non-rotatin g reference frame (LNRF) which is the orthonormal frame fixed to the 
coordinate and is called as a fixed frame in Park and (3) the comoving frame where the interactions between the matter 

and the radiation are introdu ced. As a fixed tetrad which is the orthonormal tetrad fixed with respect to the coordinates, we 
use the tetrad for the LNRF i Bardeen. Press fe Teukolskv 1972 : Frolov fc Novikov 1998h . This frame is also cal led as a zero 
angular momentum observer (ZAMO) frame. This frame corresponds the generalization of the fixed frame used in IParkl (|2006h 
for the space-time around a rotating black hole. When using the Boyer-Lindquist coordinate, the LNRF is the frame where 
the fiducial observer in this frame is rotating around the black hole with the angular velocity ui of the frame-dragging due to 
the black hole's rotation with zero radial velocity. The base e& = d/dx a can be expressed by the coordinate base d/dx 11 as 



d_ = 1 ( d_ _ g4 ,d_ \ d___ 
9* a \ at d<p)' df 



1 d 

dr ' 



d 



d 



d 



a 



(10) 



Here, (a; , x 1 , x 2 , x 3 ) = (t, f, 9, 4>). 

In the LNRF, a fiducial observer who is fixed with respect to the coordinates see the fluid with the three velocity whose 
components are calculated as 



V = -r, (i = r, 



(11) 



where u a is the four-velocity in the LNRF. The hat denotes the physical quantities measured in the LNRF. The components 
of the three velocity are explicitly calculated as 



(12) 



where f2 = jvt is the angular velocity and we have used u 
calculated as 

1= (1 



-Uf = au f . The Lorentz factor 7 for this three velocity is 



~2\-l/2 t 
V = OLU 



where v 2 = v ■ v = vtv 1 = v 2 + vj, + v 2 ,. 

A tetrad base for the comoving frame d/dx a is calculated by the Lorentz transformation as 



a 

dx° 



A^(v) 



d 
dx0'' 



(13) 



(14) 



where the bar denotes the physical quantities measured in the comoving frame. The components of the Lorentz transformation 



A a ,(v) are given as A* f = 7, A l f = jv\ A* 5 = jvj and A'^ = 5*. + v'vjj 2 /(l + 7) (i, j = r, 



Here, (a; , x 1 , x 2 , x :i ) — 



(t, r, 9, (j>). The components of the base of the comoving tetrad d/dx a can be expressed by the coordinate base d/dx a as 
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9 jv r 9 1 

df a dt 



d<j> a dt 






7tv 


a 


7+1) 






7^9 

a 




J+l) 


90' 


f \ 

7 + iy 


] 9 



(15) 



In the same way, the inverse transformation from the base of the comoving tetrad to the coordinate base is calculated by 
using the inverse Lorentz transformation A°^(— v), and the coordinate base d/dx a is calculated from the base of the comoving 
tetrad d/dx a as 



ld_ 

a dt 



9 



^jdi + 



— JV r H VrVtf, 



7 2 
7 + 1 



9 



-yve H 



7+1 



9 



+ % -7*0 + 



1 1 -2 

1 + «c* 



7+1 



1 


9 


\/jrr 


9r 


1 


9 


\flee 


90 


1 


9 



1 + &? 



9 

9 „ „ 
' 1Ve dt + VrVe 



7 2 
7 + 1 



9 „ „ 
■sz + ivve 
or 



7 2 
7 + 1 



7 2 



9_ 

df 



1 1 -2 
1 + u» 



7 + 1 y 9r 



9 „ „ 
— + «ew^> 



7 2 
7+1 

7 2 
7+1 

£.2 



9 „ „ 



9 . „ 



7 + 1 / 90 



9 

7^ + 



i 1 -2 
1 +V4 



d_ 

7+ V 90' 

_9_ 

90' 

_9_ 

90' 



7 2 



7 2 
7 + 1 



(16) 



3.2 Radiation moments 



The radiation energy density E, the radiation flux F l and the radiation pressure tensor P 1J are defined as the zeroth, the first 
and the second moments of the specific intensity I v {x tl , n), respectively. We denote the radiation moments as 



E 



iodvdtl, F i 



when measured in the LNRF, and 



IpdMQ, F 



dudn, p 13 = 



Iorth 3 dvdCl, 



(17) 



(18) 



when measured in the comoving frame. Correspondingly, the radiation stress tensors for the LNRF and the comoving frame 
are given as 



R a/3 = 



and 



E 


pr 


F e 


p<t> 


pr 


prr 


pre 


pr<l> 


F e 


pr0 


pee 


pe<t> 


p4> 


pr<t> 


pe<t> 


p<t"l> 


E 


pr 


po 


p4> 


pr 


prr 


pre 


pr<t> 


F e 


pr6 


pee 


pe<p 


p<t> 


pr<t> 


pe<t> 


p<t"l> 



(19) 



(20) 



/ 



respectively. The contravariant components of the radiation stress tensor R al3 are calculated from R al3 by the transformation 
as 



a dr a Br 13 - - 

pa/3 _ p^ v 

dxi 1 dx° 



(21) 
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and explicitly given as 



R ap = 



1 

a 



E_ 

s\'2 



OL\J% 

F 



ay lee 

p'i> 



'1H 



E 



prr 
Irr 



pr<t> 



1 

lee 



prt 



y/lrrlee 

pee 
fee 



'1U 



1 

a 



a 



jrd> 



'It 
1 



'144 
pe<t> 



lee \ V 7<W 



2^ 



, F 



1<P4> a 



F* + 



E 



By using the Lorentz transformations, the radiation moments measured in the comoving frame are calculated from those 
measured in the LNRF as 



R 



a/3 



and explicitly given as (|Mihalas fc Mihalas! Il984l ; IParklbood) 

E = 7 2 (E-2v r F + v l v J P' j ) , 



F % = -7«'E + 7 

pij 



Sj+ I 7 + — : I vv 1 



7 + 1 



F 3 - ivj S, 



I 2 
7 + 1 



I 2 



(23) 

(24) 
(25) 

(26) 



YFv'E - 7 + V51 + 2-^-jWv k j F k + \Sl + Y^vv k j + "r^y^J P ■ 

Inversely, the radiation moments measured in the LNRF are calculated from those measured in the comoving frame as 
R && =A%(v)A p (v)R* D , (27) 
and explicitly given as 

E = -f (E + 2v i F i +v i VjP ii ) , (28) 



F 1 = 7i'£ + 7 



s' + h + 



7 + 1 



F' + V, i^ + ^Y^fe)^, 



7 Vfi J P + 7 ( -C"^ + &51 + ]F k +\^5l + Y^vv k j[Si + ^TJ^^ j P 



I 2 



7 + 



I 2 



7 + 



(29) 
(30) 



3.3 Radiation four-force density 



The radiation four-force density measured in the comoving frame is given as (Mi halas fe Mih alas 1984) 



G a = - 

c 



(xF — F)n a dvdQ., 



(31) 



where \ an d fj are the mean opacity and the emissivity measured in the comoving frame, respectively. The time component 
G 4 has the dimension F 1 times the net rate of the radiation energy per unit volume, and the spatial component G l has the 
dimension of the net rate of the momentum exchange between the matter and t he radiation. The time component o f the 
radiation four-force density measured in the comoving frame can be calculated as jMihalas fc Mihalas! [liii ; IParklbOQrj ) 



G = T - A 

where the heating function T and the cooling function A are defined as 



XlpdvdQ, A 



rjdvdQ,. 



The components of the radiation force G a is calculated from those in the comoving frame G a by the transformation 

a _ dx a a 



(32) 



(33) 



(34) 
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and explicitly given as 



G' 

G 6 
G'" 



7/ 
a V 



G + v t G l , 



1 



7rr- 
1 



G f +JV r G f + ^— -VrViG* ) , 

7 + 1 



7 2 



G° + yveG* + ^l—veViG* 
lee \ 7 + 1 

-^L+7^-^^ + 7 
\/7<«> 



7«<A 



'7w(l + 1) Q 



mG\ 



(35) 



3.4 Radiation hydrodynamic equations 

3.4-1 Continuity equation 

As a continuity equation, now we consider the particle number conservation, (nu OL )- a — 0. This equation can be calculated as 
d ( t r-\ d ( t r \ , (36) 



where 7 = det7sj. In the case of the Kerr metric written by the Boyer-Lindquist coordinate 7 = (EA/A) sin 2 8. 



3.4-2 Hydrodynamic equations 

The relativistic Euler equations are obtained by the projection of the equation of the energy-momentum tensor = G a 

on the specific directions by using the projection tensor P a ^ — g a/3 + u a u' 3 as PpT l3S s = PpG 13 . From this, we can obtain 
p o hgu a . u d + g a(, P g ,p + u a u P g , f 3 =G a + ^upG 13 . 

The radial part of the Euler equation, i.e. momentum conservation for r-direction, is calculated as 



t du r , i du r Poh R 
Pohg u—+ Po h s u— + — 



(ln7rr) iT (u r ) 2 + 2 Qn^rr)^ M ' u ~ ( m 7ee), r ( u Y 



[ w 0n 7 *) jP (fi + + 2 W( 3? r (n + /3 ) - 2a 2 (Ina) >r ] + 



= -a'tfvTG 1 + \l + irr(u) 2 } G r + ~ieeuu 9 G e + 7w (fl + ,S*)(G* + ^G^JtiV, 



+ u I it 



9* 



+ u 



(37) 



where the derivatives for In a, /?^, ln 7rr , In "/ee and In 7< ^ with respect to r and # are given in App.|X] The terms including /3*'(= 
— u)) explicitly show the effects of the frame-dragging due to the black hole's rotation. Here, we have used the relation gu t r(u t ) 2 + 
2 5t0 , r itV + g^A^f = (V) 2 [7^(ln7^), r (^ + /3' / ') 2 + 27^At(" + ^)-2a 2 (lna) ir ]. This terms are related to the 
centrifugal force. Actually, this relation can be also calculated as gu,r( ut ) 2 + 2g t )r u t u' t ' + g^0,r(w^) 2 = (ln^^rT^^U 1 ) 2 ^ — 
n_)(fi — fi+) where fi± is calculated as 



9t<t>,r 

g§4>,r 



9t<f>,: 



9tt,; 



±y/m/ sin 6 



rS 2 
2r 2 - : 



(38) 



m(r 2 + a 2 — E) + ma 2 sin 2 6 ± y^a sin 2 



These angular velocity Q± is the generalization of the Keplerian angular velocity = ±m 1/ ' 2 /(r 3 ' /2 ± am 1 ' 2 ), and in the 
limit of 6 = 7r/2, we obtain fi± = fi^. In the right hand side of Eq. (I37|l . the terms not including the derivatives represents 
the gravitational acceleration and the centrifugal acceleration, respectively. 

The Euler equation in ^-direction, i.e. momentum conservation in ^-direction is calculated as 



, t du , idu pohe 

poK u — + P0 h,u — + — 



p h s (1 



(111799) . (u) 2 + 2 (In -fee) uu - (ln 7rr ) — {u) 2 



[7** 0*7*),, (0 + /? ) 2 + 27^^(0 + /3*) - 2a 2 (Ina),,] + + u 9 (V^f + 



796 



-aVu'G* + t^mVCT + 1 + 7ee (^) 2 G H + ^4,(0, + j3 <p )(G' p + ff"G,)u u 



(39) 
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In the same way, the Euler equation in (^-direction, i.e. momentum conservation in 0-direction is calculated as 



pahgU + pohgU -p—r + pohgU u U}l + [V) 
+ poh g u u <^(fi + 



(In 7^) a + fit 



a 



2/3* (Ina) iS +j8f s } + 



-2/9* 0no) iP +/3*} 

i dP 3 , /3* / 9P 9 



+ 1 — I - — - - p 

7<w dcj> a 2 V 9t 



^P 9 



+ u*^) = -a 2 uVG* + 7 r^ r u^G r + 7ee uVG 9 + G* + 7w (fl + /3*)(G* + /3*G*)uV 



(40) 



The for the local energy conservation is obtained from u a T a ^j = u a G a calculated as 
dt \ n J dx l \ t 

The right hand side of Eq. (|41[) is also calculated by the heating and cooling function defined in the comoving frame as 
u a G a = -G l = A — F. (42) 



,„ t -_*_ f £^^ +u t^ +u >^JL = _ q VG< + 7 ^ r G r + 7efl ^G 9 +7^(fi + /9*)(G* + p*G l )u t . (41) 



3.4.3 Radiation moment equations 

The radiation moment equation il = — G a gives the equation for the energy density, the radiation flux and the radiation 
pressure tensor. The radiation energy equation is obtained from f-component of this equation R ta . a = — G* calculated as 



eft a-^/7 



+ -rr + (ff gtt, r + g v gtj>,r) R + (g gt*,r + g 



t<i> 



R r 



( tt . t<4 \ r>t9 , I tt , tib 

+ (g gu,e + g v gt<t,,e)R + (g gt4,.e + g^g 



)R U<P = -G\ 

The radiation momentum equation in r-direction P r< * a = — G r is calculated as 



(43) 



dR tr 1 
+ 



9t 0^7 
1 



d_ 

dr 



97? * 1 rr / rr r 9 g 

H ^: 1" TTff (ffrr.r-H +2g rr . e R -gee.rR 

06 2 v 



<T (Stt.r-R" + 2 fl ^,,.P t * + g^,rR^) = -G 



(44) 



In the similar manner, the radiation momentum equation in ^-direction 7? ™ a = — G e is calculated as 

dP e * 1 9fl . 



-\g BB (gtt,eR u + 2g tct> ,eR t4 ' + g^fiR 94 ') = -G' 



d<t> 



■)R 0e + 2gee, r R rB - g r r,eR rr ) 



(45) 



The radiation momentum equation in (^-direction R 4a , a = — G* is calculated as 



dR u 
~dt 



+ 



(HI 



dR <t>4> 



+ —k— + (g^gt^.r + g^gu.r) R tr + (g^g^.r + g^g^r) R r 



+ (g 9tt,e + g gtt.e) R + {g g<t>4>,e + g gt^.e) R = -G 



(46) 



Finally, by inserting the components of the radiation stress tensor given by Eq. (|22|l in t o Eqs. (1431) . (144 |l . (|45|l and (|46|l . we 
obtain the radiation moment equations corresponding to Eqs. (49), (52), (53) and (54) in I Par id (|2006l ) as 



d_ f E_\ 1 d 

dt \a 2 J ay^dr \]J j rr 

+2(lna) r - +2 (In a 



-l-F r ) + 1 9 



7 98 \ \J jeo 



a 



■pe- 
at 



P% 1 144 pr<t> 1 144 pH _ Qt 

a J -yee q2 V > r ^ \ l ef > 



d_ 

dt 



F r 



a ylr 
On a). 



£ ■ 



1 d ( r P r ' 



F" t> +~(kij rr ) r - — + (ln 7rr ) 



a V 1 





19/ r- P r ® 
a V 7- 



«V7 



9^ 



a 



'IuPTt 



"~Yrr 



d_ 

dt 



Offrr 

1 d 



\J lrr"/ee J \flrr \ \flH, 

P rt> Qniee),rp 9e _ (In 7*0 ), r = _ G ,. 



' 7 rr 
pre 



\J Irrlee 2 >'' 



2 7r 



cty"yee I «VT 



9r 



1 a / r P Be \ d 
, — = 1 + — oy 7 — 1 + — 

y/lrrlee J aiy/jOV \ lee 



lee \ \/l<i><t> 



a 



(47) 



(48) 
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(In a) „ „ 
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6 £,</> 



a^ee 






prr 

9 lee 


1 d 


a \fl 


ty/l dr 
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+ d<j> 






+ ^ (ln7ee), 9 5 (In 7<M>), 9 — 



-G . (49) 



7<W 



+ [/# - 2j8*(lna), P ] — ^= + [/$ - 2/J*(Ina),«] — ^= + 
oAr ' oid iee 



(In- 



+ 



(In- 



pS</> 



-G , 



+ 



V Irrlqxf, 



\fieelH, 

where the derivatives for In a, /3 , ln7 rr , In lee and In 7^ with respect to r and 9 are given in App. 1X1 



(50) 



4 CONCLUDING REMARKS 

As is well known, the moment equations which are truncated at the finite order do not constitute a complete system of 
equations. That is, the number of equations are s maller than the number of variables to be solved. Thus, we need the additional 
equa t ions to close the system of the equ ations 1 Chandrasekhar 1960l ; Mihalasl 197p| ; Pomranine 19731 ; Rvbicki fc Lightman 
19791 : iMihalas fc Mihalasl [l98i : Ishulll99lh such as the Eddington factors. In the diffusion limit where the photon mean-free 
path is much smaller than the characteristic length of the system, i.e. in the optically thick region, the so-called Eddington 
approximation P %3 = (<5 lJ /3)E is valid. On the other hand, in the optically thin region, the radiation becomes anisotropic in 
general. So, in such case, the full angle-dependent radiative transfer equations should be solved, and sometimes the variable 
Eddington factors are introduced. Some authors have expressed the variable Eddington factors as a function of the optical 
depth. Additionally, it is pointed out that the influence of the flow velocity should be also taken into account w hen constructing 
the c l osure relations. These prob l ems and their calculation methods are extensively studied in the past studies |Auer fc Mihalasl 
1970l: Irlummer fc Rvbickilll97ll: fTamazawa et al.lll975l: IMihalas! fl98ol; IMihalas fc Mihalasl 1 1984: Igchinder fc Bludmanlll98S 



Nobili. Turolla fc Zampierilll993l : lYin fc Millerl Il995l : iKato, Fukue fc Mineshigelll998l : iFukuef 20061 ). The similar analysis for 
the closure relations will be studied in the case of the Kerr space-time in future. In the present study, we use the Boyer- 
Lindquist coordinate which is a frequently used coordinate in the astrophysics around a rotating black hole. However, there 
is a coordinate singularity at the horizon in this coordinate. Future studies also resolve this singularity by using the another 
coordinate with no singularity at the horizon (see, e.g. Takahashi 2007). 

In this study, we have derived the radiation hydrodynamic equation in the Kerr space-time. While the interactions 
between the matter and the radiation are defined and calculated in the comoving frame, the derivatives which are used to 
describe the global evolutions of both the matter and the radiation are calculated in the BLF. Both the matter and the 
radiation are influenced by the frame-dragging due to the black hole's rotation. In this approach, as a fixed orthonormal 
frame, we use the so-called locally non-rotating reference frame (LNRF) which is firstly calculated in the Kerr space-time by 
Bardeen. Press fc Teukolskv ( 19721 ) by using the Boyer-Lindquist coordinate. The special relativistic effects such as beaming 
effects are introduced when the Lorentz transformation between the LNRF and the comoving frame. On the other hand, 
the general relativistic effects such as frame-dragging effects due the black hole's rotation and the gravitational redshift are 
introduced by the tetrads describing the BLF. 
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APPENDIX A: METRIC COMPONENTS AND DIFFERENTIAL VALUES 



Here, we present the explicit expressions for the metric components and the their differential values with respect to r and 
used in this paper. Non-zero components of the metric are given as 



gtt = -a + fop? = - 1 1 — 



g rr = j rr — — , g ee = j 0g = E, 



and 
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'ea' 



tip 
9 = — 



gt<j> = j3<f, = 



2mar 



2mar sin 9 



-, 944 = 7^ = 



A sin 2 9 
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; 7 



ee 

9 =7 



ct 

ee 1 
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<i>\2 



A sin 



m 2 e i 1 



2mr 



A 

E ' J ' E' 
The differential values of the metric used in this paper are given as 



(lna) r = i [(lnE), r + (lnA) r - (lnA) „] , 0* (^A) r 

(ln 7rr ) ir = (lnE) :r + (lnA) r , (ln 7ee ) r = (In E) r , (In 7^ 
(lna) i9 = \ [(lnS) i9 + (In A) , - (lnA),,] , = w (In A\„ , 



(lnA) r -(lnE) rJ 



(ln7rr), s = (lnE) <e + (InA), e , (In-yeo)^ = (lnE)_ e , (ln 7 ^) i9 = (lnA) e + 



tan 9 



(InE), 



where 



(lnS) iP = J, (lnA) <r = 2(r A m) , (lnA) ir = |[(r-m)E + (r + m)(r 2 +a 2 )], 



(A: 



(A: 



(A: 
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2 A 

(In T,) e = —a 2 sin 20, (lnA) j9 =0, (In A) e = sin 20. 

Here, we also give the differential values of the metric as 



(A4) 



gu,r 



and 



gtt,e ~- 
grr,e 



2m 
~E~ 

2 
A 



1 - 



2ma sin 2 / 2r : 



E 

(r-m)S 



9t<p,r 



, gee,r = 2r, 



~ 1 ) > ff^.r = 2sin 



r — m + 



(3r 2 +a 2 )m 2mr 2 (r 2 + a 2 ) 
S 



(A5) 



2ma 2 r . „„ 2mar(r 2 +o 2 ) . „„ 

E2 sin 29, 3*0,9 = — sin 20, g^e 

a 2 

sin 20, = —a 2 sin 20. 



A + 



2mr(r 2 + a 2 ) 2 
S2 



sin 20, 
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